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We present a measurement of the ratio of top-quark branching fractions R = Bit — > Wb) /B(t — > 
Wq), where q can be a b, s or a d quark, using lepton-plus-jets and dilepton data sets with inte- 
grated luminosity of ~162 pb _1 collected with the Collider Detector at Fermilab during Run II of 
the Tevatron. The measurement is derived from the relative numbers of ti events with different 
multiplicity of identified secondary vertices. We set a lower limit of R > 0.61 at 95% confidence 
level. 

PACS numbers: 14.65.Ha, 12.15.Hh 



The top quark as described by the Standard Model 
(SM) is expected to decay to a W boson and a bot- 
tom quark at least 99.8% of the time at 90% confidence 
level (CL) Q. The Cabbibo-Kobayashi-Maskawa (CKM) 
quark- mixing matrix 0,0 element |Vtb| is expected to 
be very close to unity from the assumption of a unitary, 
three-generation matrix and the measured small values 
of \Vub\ and \V c b\ 0- A measurement of the ratio of top- 
quark branching fractions R = B(t — > Wb)/B(t — > Wq), 
where q can be a b, s or a d quark, significantly less than 



unity would contradict our current theoretical assump- 
tions, implying either non-SM top decay, a non-SM back- 
ground to top-pair production, or a fourth generation of 
quarks. A previous measurement has set a lower limit of 
R > 0.56 at 95% CL 0. In this Letter we present a mea- 
surement of R using ti events collected at the Collider 
Detector at Fermilab (CDF) during Run II of the Teva- 
tron, a proton- antiproton collider with center of mass 
energy of = 1.96 TeV. The integrated luminosity of 
the data sample used in this analysis is ~162 pb _1 . 
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Our measurement uses tt-pair events. The lifetime of 
top is too short for hadronization to occur, and the SM 
strongly favors an essentially immediate decay of each 
quark to a real W boson and weak-isospin -1/2 quark; if 
R = 1, this is always a b quark. To maintain high detec- 
tion and trigger efficiencies and low background levels, 
we only consider ti final states in which at least one W 
has decayed leptonically. Events in which one W decays 
lcptonically are called "lepton-plus-jets" (L+J) events, 
and events with two leptonic decays are called "dilepton" 
(DIL) events. Values of R are determined separately for 
each of these sets of events, and are combined in the end 
to set a lower limit on R. The greater statistical power 
comes from the L+J sample. 

The measurement requires both the counting of b- 
quark jets and the determination of the ti content as a 
function of the 6-quark multiplicity. We identify ( "tag" ) 
6-quark jets by identifying displaced secondary vertices 
using the SECVTX algorithm R is extracted from 
the relative rates of events with zero, one, and two tags; 
any two rates determine R uniquely, while all three rates 
jointly overdetermine R. A novel feature of this mea- 
surement is the inclusion of the 0-tag L+J event rate, 
which is determined using event kinematics and an ar- 
tificial neural net (ANN) technique. As R depends only 
on relative rates, this measurement is independent of any 
assumptions of the overall ti cross section. However, our 
measurement of R does depend critically on the knowl- 
edge of the efficiency to identify b jets. To extract R we 
use the efficiency to tag jets in ti events estimated with 
a Monte Carlo (MC) sample in which tagging efficiencies 
have been tuned to match jet data 0- 

The CDF detector for Run II [f| consists of a charged- 
particle tracking system in a magnetic field of 1.4 T, seg- 
mented electromagnetic and hadronic calorimeters and 
muon detectors. A silicon microstrip detector provides 
tracking over the radial range 1.5 to 28 cm, and is essen- 
tial for the detection of displaced secondary vertices. The 
fiducial region of the silicon detector covers the pseudo- 
rapidity range \r)\ < 2, while the central tracking system 
and muon chambers provide coverage for \r)\ < 1 |2|. A 
three-level trigger system is used to select events with 
electron (muon) candidates with Et (jpt) > 18 GcV (18 
GeV/c), which form the data set for this analysis. 

L+J events consist of one isolated high-pr lepton (e or 
fi), large missing transverse energy ($ T ) due to the un- 
detected neutrino, and four hadronic jets. Two of these 
jets arise from the hadronic decay of the other W, and 
the other two arise from the top-daughter quarks q. The 
L+J selection requirements are described in detail else- 
where 0- Briefly, we require the presence of an isolated 
lepton which has transverse momentum greater than 
20 GeV/c, that $J T is at least 20 GeV, and that there is a 
minimum of four jets, clustered with a cone-based algo- 
rithm having cone dimension AR = y/ (A0) 2 + (A77) 2 = 
0.4, within |t;| < 2 and with corrected transverse en- 



TABLE I: Summary of observed number of events with 
i tags in the L+J and DIL samples, with estimates of nominal 
ti event-tagging efficiencies, background levels and expected 
event yields. The L+J 0-tag background is measured with an 
ANN. The efficiency estimates and the 1-tag and 2-tag L+J 
background estimates are given for R = 1. Equations and 
121 are used for the calculation of the expected total number 
of events N^ p . The statistical and systematic uncertainties 
have been combined. 



Lepton + Jets (L+J) 


O-tag 


1-tag 


2- 


■tag 


Efficiency (<h(R= 1)) 


0.45 ± 0.03 


0.43 ± 0.02 


0.12 


± 0.02 


Background (A^ kg ) 


62.4 + 9.0 


4.2 ± 0.7 


0.2 


±0.1 


Total expected (N-* p ) 


80.4 + 5.2 


21.5 + 4.1 


5.0 


± 1.4 


Observed (N° bs ) 


79 


23 
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Dileptons (DIL) 


O-tag 


1-tag 


2- 


■tag 


Efficiency (e 4 (it= 1)) 


0.47 ± 0.03 


0.43 ± 0.02 


0.10 


± 0.02 


Background (A^ kg ) 


2.0 + 0.6 


0.2 + 0.1 


< 


0.01 


Total expected (N-* p ) 


6.1 + 0.4 


4.0 ± 0.2 


0.9 


±0.2 


Observed (N° hs ) 


5 


4 
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ergy |5j greater than 15 GeV. These requirements select 
107 events. 

DIL events consist of two charged leptons (ee, [i\x or 
efj,), large $ T due to the undetected neutrinos, and two 
jets from the top-daughter quarks q. The DIL selection 
requirements are described in detail elsewhere Q ■ Com- 
pared to the L±J selection, we demand an additional 
lepton, but only a minimum of two energy-corrected || 
jets, with the same requirements as before. These re- 
quirements select 11 events. 

Both event samples are subdivided on the basis of the 
number of identified b jets in the event. The number of 
events in each subsample with i tagged jets are given in 
TablcQ] The 2-tag subsample is defined to include events 
with > 2 tags; in this data sample we observe no events 
with more than two tagged jets. 

In the L+J sample, the dominant background is W 
production in association with jets from QCD processes 
("W+jets" events). In the 1-tag and 2-tag subsamplcs 
we make an a priori estimate of the backgrounds with a 
collection of data-driven and simulation techniques that 
arc described in detail elsewhere [f| The backgrounds 
in these subsamplcs include W production in association 
with heavy-flavor jets (Wbb, Wcc, Wc), W production 
in association with light-flavor jets that are incorrectly 
identified as b jets ("mistags"), QCD multi-jet ("QCD") 
events containing fake or real leptons and/or incorrectly- 
measured $ T , dibosons (WW, WZ) and single-top quark 
production. The background estimate requires a small 
correction for R 7^ 1 . The background estimate for R = 1 
in these subsamples is given in Tablc|U The uncertainties 
on the estimate are dominated by uncertainties in the 
fraction of VL+jets events that include heavy flavor and 
on the normalization of the QCD background rate. 

By construction, the a priori method cannot pre- 
dict the background level in the 0-tag L+J sample, 



5 



to r 

-4— ' 

c 

CD 

a>io- 



CD 
_Q 

E 

13 



CDF Run II W+jets, 53.4 ± 7.6 

jLdt = 1 62 pb" 1 QCD multi-jet, 9.0 

■ tt, 16.6 ±7.6 

— Sum 

▼ Data (0-tag), 79 



T T T 



T T 



TO 



Q 







T T 

mm 

0.2 0.4 0.6 0.8 

Neural network output 



FIG. 1: Fit of the ANN output in the 0-tag L+J data set (tri- 
angles) with a sum of 3 components: W+jets (upper), QCD 
multi-jet (middle), and ti (lower). The QCD normalization is 
independently estimated and not varied in the fit; its shape 
is determined from the non-isolated lepton data. 



where the W+jets production rate dominates that for 
tt pairs; instead we make use of event kinematics Q. 
The artificial neural net E3 is trained with the tt sig- 
nal (HERWIG Ej) and W+jets background (HER- 
WIG+ALPGEN |l3) events simulated with a detailed 
detector description based on GEANT [lflj . There is 
an additional QCD background which is modeled us- 
ing data with non-isolated lcptons. We find optimal sig- 
nal to background discrimination with an ANN structure 
of nine input variables, one intermediate layer with ten 
nodes, and one output unit. The variables used are the 
transverse energies of the four leading jets, the minimum 
di-jet mass, the di-jet transverse mass with value closest 
to the mass of the W, the scalar sum of the transverse 
energies of all leptons and jets, the total longitudinal mo- 
mentum divided by the total transverse momentum, and 
the event aplanarity. 

The ANN output ranges from zero for background-like 
events to one for signal-like events. We perform a binned 
maximum likelihood fit of the ANN output distribution 
for the tt fraction in the 0-tag subsample. The fraction of 
events from QCD backgrounds is fixed to 11.4% in this 
fit. These events are characterized by the non-isolation 
of the lepton and small ~$ T , and the fixed rate is based on 
comparing to control regions with either low f T or poor 
isolation [5j . The resulting measurement of background 
rates in the 0-tag L+J subsample is given in Tabled The 
fit of the distribution of ANN outputs for this subsample 
is shown in Figure ^ 



Systematic uncertainties in the ANN-determined back- 
grounds are dominated by our understanding of the jet 
energy scale, the renormalization and factorization scale, 
and the shape of the QCD template and are strongly anti- 
correlated between the tt and W+jets measurements. 
Our ANN-measured tt content in the L+J sample with- 
out any tagging requirement is consistent with that found 
in our earlier measurement of the tt production cross sec- 
tion 01- T ne procedure is repeated in the 1-tag and 
2-tag samples, yielding background rates of 5.8 ± 5.2 
and O.l^o 5 respectively, consistent with the a priori es- 
timates shown in Table HJ As the a priori estimates have 
smaller uncertainties in the 1-tag and 2-tag subsamplcs, 
the ANN-dctcrmincd background level is used only for 
the L+J 0-tag subsample. 

The main backgrounds in the DIL sample arc Drcll- 
Yan production including lepton pairs from the Z res- 
onance, dibosons, and W+jets events with fake lcptons. 
The total background level of 2.2 ± 0.6 events in the DIL 
sample has been estimated elsewhere ||. The Drell-Yan 
rate in ee and events is estimated using simulated data 
normalized to the observed rate of Z events in the data. 
Other electroweak backgrounds are estimated from MC 
simulations. The fake-lepton background is estimated by 
multiplying each jet in W plus three or more jet events 
by a lepton fake rate, measured in complementary jet 
samples. 

Most of the jets in the DIL background events arise 
from generic QCD radiation. To determine the back- 
ground distribution across the i-tag subsamples, we ap- 
ply a parameterization of the probability to tag a generic 
QCD jet 0, derived from jet-triggered data samples, to 
the jets in the DIL sample, correcting for the enriched 
ti content of the sample. The resulting estimates are 
given in Table [U the background in the 2-tag subsample 
is negligible. 

The tt event-tagging efficiency e^, defined as the proba- 
bility to observe i tags in a ti event, depends on the fidu- 
cial acceptances for jets that can potentially be tagged, 
and the efficiencies to tag those jets . Those efficiencies 
in turn depend on the species of the underlying quark in 
the jet. The efficiency e, ; depends strongly on R, as R ^ 1 
implies fewer b jets available for tagging, and more light- 
quark jets available instead. We use the jet acceptances 
and tagging efficiencies to parameterize ti(R). These 
quantities are estimated with a sample of simulated ti 
events from the PYTHIA [l]| generator and CDF de- 
tector simulation, and their uncertainties are dominated 
by our understanding of the control samples of jet data 
used to calibrate tagging efficiencies in the simulation. 
The leading determiner of is the efficiency to tag a b 
jet from the decay t — » Wb\ e b = 0.44 ± 0.04 for b jets 
falling within the fiducial acceptance and having at least 
two tracks with silicon information. The e, values also 
have small contributions from the efficiencies to tag jets 
from W — > cs hadronic decays and from additional QCD 
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radiation in tt events. The nominal values of ej for R = 1 
are given m TablcUl The value of e (e 2 ) changes by -0.28 
(0.09) as R changes from 0.5 to 1. 

The expected event yield in each of the three tagged 
subsets of each of the L+ J and DIL samples is 



Nr = K c -e,(R) + N„ 



(1) 



where 7V bkg is the number of background events in the 
i-tag subsample and N-£ c is an estimate of the inclusive 
number of ti events in the sample, determined by 



N: 



N 



bkgx 



(2) 



where iV° bs is the observed number of events in each sub- 
sample. In this construction, the measured value of R 
is independent of any assumption of the overall rate of 
ti production, and is thus sensitive only to the relative 
numbers of ti events with i tags. 

The full likelihood is a product of independent likeli- 
hoods for the L+J and DIL samples. Each likelihood is a 
product of Poisson functions comparing A° bs to N^ xp for 
each value of i, multiplied by Gaussian functions which 
incorporate systematic uncertainties in the event-tagging 
efficiencies and backgrounds, taking into account the cor- 
relations across the different subsamplcs. These include 
correlations in the event-tagging efficiencies through the 
single-jet tagging efficiencies; in the common methodol- 
ogy of the a priori estimates in the tagged L+J samples; 
and in the overall normalization of the DIL backgrounds. 
There are a total of five free parameters in the likelihood 
to account for these systematic uncertainties. 

The resulting likelihood as a function of R is 
shown in Figure |2 along with the negative loga- 
rithm of the likelihood. We find a central value of 
R = 1.12tp;ig(stat)t°; 1 3(syst). The dominant system- 
atic uncertainties arise from the uncertainty on the back- 
ground measurement in the 0-tag L+J sample (jZo^ii) and 
from the overall normalization of the tagging efficiencies 
(looc!)- Taken separately, the two final states of tt give 
consistent results for R; the L+J sample alone yields 
R = 1.02;to 2o±o.i3> and thc DIL sample alone yields 
R = 1.411°' loio i3- These R results are consistent with 
the SM expectations. 

The ratio R can only take on physical values between 
zero and unity. We use the Feldman-Cousins prescrip- 
tion 0| to set a lower limit on R. We generate ensem- 
bles of pseudo-experiments for different input values of 
R (Rtrue), and vary the input quantities of the analysis, 
e.g. the background estimates, taking correlations into 
account. Using the likelihood-ratio ordering principle, we 
find the acceptance intervals as shown in Figure [21 With 
our measured value of R, we find that R > 0.61 at the 
95% CL. 

Our lower limit on R is the strongest limit on this 
top-quark branching ratio to date. Within the SM, R = 
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FIG. 2: The upper plot shows the likelihood as a function 
of 7? (inset) and its negative logarithm. The intersections of 
the horizontal line ln(L) = —0.5 with the likelihood define 
the statistical la errors on R. The lower plot shows 95% 
(outer), 90% (central), and 68% (inner) CL bands for R tr ue 
as a function of R. Our measurement of R = 1.12 (vertical 
line) implies R > 0.61 at the 95% CL (horizontal line). 



\v tb 



Vtb ^+lVt P+MdP ' U P ^° phase-space factors. Assuming 
three generations and the unitarity of the CKM matrix, 
the denominator is unity, and we estimate \Vtb\ > 0.78 
at 95% CL. All of our measurements of R are consistent 
with thc SM expectations. 
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